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Abstract 

In this paper, we present a new approach to the discrete version of the Wormlike Chain Model 
(WCM) of semiflexible polymers. Our solution to the model is based on a new computational 
technique called the Generalized Borel Transform (GBT) which we use to study the statistical 
mechanics of semiflexible polymer chains. Specifically, we evaluate the characteristic function of 
the model approximately. Afterward, we compute the polymer propagator of the model using 
the GBT and find an expression valid for polymers with any number of segments and values of the 
semiflexibility parameter. This expression captures the limits of flexible and infinitely stiff polymers 
exactly. In between, a smooth and approximate crossover behavior is predicted. Another property 
of our propagator is that it fulfills the condition of finite extensibility of the polymer chain. We 
have also calculated the single chain structure factor. This property is a decreasing function of 
the wave vector, k, until a plateau is reached. Our computations clearly show that the structure 
factor decreases faster with increasing wave vector when the semiflexibility parameter is increased. 
Furthermore, when the wave vector is large enough, there is a regime where the structure factor 
follows an approximate power law of the form k~ e even for short polymer chains. is equal to two 
for flexible polymers and to one for rigid chains. We also compare our results to the predictions of 
other models. 
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I. INTRODUCTION 



In recent years, experimental studies of biological macromolecules have motivated intense 
research in the field of statistical mechanics of single semiflexible polymer chains. Indeed, 
studies like force-elongation measurements of different biological (DNA, titin, tenascin) and 
synthetic (polyethylene glycol, polyvinyl alcohol) polymers using Atomic Force Microscopes, 
Optical Tweezers and other recently developed tools are abundant The measured force- 
elongation curve is generally fitted to the prediction of the Wormlike Chain Model (WCM) of 
semiflexible chains^, originally proposed by Kratky and PorodQ]. From this fit, parameters 
like the persistence length of the biopolymer are extracted. Another kind of experiments has 
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. These properties are determined 



targeted the mechanical properties of eukaryotic cells 
by an assembly of protein fibers called the cytoskeleton. This three dimensional assembly is 
made of the cytoskeletal polymers (microtubules, actin filaments, etc.). All these polymers 
are semiflexible polymers at the relevant length scales (a few microns at most). Thus, the 
predictions of the WCM are very relevant for the understanding of the physical behavior of 
cytoskeletal polymers. 

The WCM was originally proposed by Kratky and Porod in 1949^1 and reformulated 
using field theoretic methods by Saito et al\^ in 1967. In this model, the polymer chain 
displays resistance to bending deformations. This resistance is modeled using a free energy 
that penalizes bending the polymer backbone. The free energy depends on parameters (elas- 
tic constants) that are a consequence of many short-range monomer-monomer interactions. 
Explicitly, the free energy is 

' 'd 2 R(s) N2 



where R (s) is the vectorial field that represents the polymer chain, s is the arc of length 
parameter, L is the contour length of the polymer and k is the bending modulus. In addition, 
the local inextensibility constraint |dR(s) / ds| = 1 must be satisfied. As a consequence of 
the bending rigidity , a wormlike chain is characterized by a persistence length (proportional 
to the bending modulus) such that, if the length scale is shorter than the persistence length, 
then the chain behaves like a rod while, if the length scale is larger than the persistence 
length, then the chain is governed by the configurational entropy that favors the random- 



walk conformations. 

The local inextensibility constraint has not allowed researchers to find an exact solution 
to the WCM. Indeed, the constraint |dR(s) / ds| = 1 is written using a Dirac delta distri- 
bution in infinite dimensions. Depending on how the constraint is written, |dR(s) /ds| = 1 
or (dR(s) /ds) 2 = 1, we get an Edwards Hamiltonian that is non- analytic or non-linear, 
respectively. Consequently, there is no exact solution of this model at present. However, a 
few properties like the first few moments of the distribution of the end-to-end distance^, 0| 
are known exactly. 

The aforementioned complexity of the WCM has motivated many approximate treat- 
ments of semiflexible polymers. For example, Fixman and KovacQ developed a modified 
Gaussian model for stiff polymer chains under an external field (external force). In this 
approach, they computed an approximate distribution for the bond vectors from which they 
were able to compute the partition function and average end-to-end vector. An alterna- 
tive approach was proposed by Harris and Hearst [S] who developed a distribution for the 
continuous model from which they were able to compute the two-point correlation function 
and, consequently, the mean-square end-to-end distance and radius of gyration. Another 
statistical property of interest has been the distribution function of the end-to-end distance 
or its Fourier transform. Many approximations for this function have been proposed. For 
example, different expansions of the distribution in inverse powers of the number of segments 
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have been developed[9, lOj. Similarly, perturbations with respect to the rodlike limit have 
been derived [111. Other approaches to the distribution function of the end-to-end distance 
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of semiflexible polymers have led to modified Gaussian functions [12, LL2fl. Finally, many 
approximations have been proposed for the structure factor Q|. 

As stated in the previous paragraph, most of the approximate treatments of the distri- 
bution function of the end-to-end distance have been perturbative in nature. Indeed, these 
approximations were perturbation expansions with respect to the flexible or rigid chain lim- 
its. A different approach to semiflexible polymers was taken by Kholodenko QHQ. In 
this model, the Euclidean version of the Dirac propagator is used to predict the conforma- 
tional properties of semiflexible polymers. In particular .the single chain structure factor has 
been used to describe experimental data quantitatively (l^fl. Recently, Winkler has proposed 
another treatment of semiflexible polymers j^. In this work, an approximate expression for 



the distribution function valid for any value of the stiffness of the polymer backbone was 
developed using the Maximum Entropy Principle. 

Motivated by the experimental studies done on semiflexible polymers and our incomplete 
understanding of the properties of the WCM, we have developed a new approach that 
captures many physical properties of the model, like the limits of flexible and rigid polymers, 
exactly and provides approximate crossover behaviors for all the distribution functions. To 
accomplish this goal, we have employed a computational technique called the Generalized 
Borel Transform (GBT) which was taken from Quantum Mechanics and Quantum Field 
Theories (21, [22, 23|. This method computes Mellin/Laplace transforms exactly. We provide 
a brief summary of this technique in Appendix A. 

This paper is organized as follows. In section II, we evaluate the characteristic func- 
tion of the WCM approximately such that some physical constraints are satisfied exactly. 
Furthermore, we evaluate the distribution function (polymer propagator) using the GBT 
and compute the single chain structure factor. In Section III we discuss the results of our 
calculations which are valid for any value of the semiflexibility of the polymer. Section IV 
contains the conclusions of our work and some speculations about applications to polymer 
physics of possible extensions of the GBT. The details of the mathematical calculations are 
presented in the appendices. 

II. THEORY 

A. The Model and Evaluation of the Characteristic Function 

Consider a polymer chain modeled as a sequence of n bond vectors (Ri, R 2 , R n ) con- 
nected in a sequential manner. In addition, let us assume that the length of each bond 
vector is / (=Kuhn length) and that pairs of consecutive bond vectors try to be parallel to 



each other. This pre 
following expression 



'erential orientation is modeled with a Boltzmann weight given by the 
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where a is the strength of the interaction in units of thermal energy (= kn T). Inserting Eq. 
(J2J into the expression for the propagator of the Random Flight Model[25], we obtain the 
following expression for the polymer propagator of semiflexible chains 

Tl f TV \ / Tl — 1 \ 

P (R, n,a) = j d {R k } f[ r (R,) S f £ R i ~ R j ex P [® E R +i ' R * J > ( 3 ) 
where R is the end-to-end vector and r (Rj) is given by the formula 

The propagator, Eq. (JHJ), is not normalized. 

We proceed to express the delta function using its Fourier representation [3| then, Eq. 
(JHJ) becomes 
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which can be used to define the characteristic function, K (k, a, n, I), as follows 
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3 exp (— iR • k) K (k, a, n, I). (6) 
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The mathematical expression of the characteristic function is 



j . n / n n— 1 \ 

tf(k,a,n,Z) = — — s / d{R fc }n^(|R-i|-Oexp U'E R i- k + «E R wRi I ^ 

I 47 ™ J J j=l \ 3=1 i=l / 

Note that since the polymer propagator, Eq. Q, is not normalized then the characteristic 
function, Eq. 0, does not approach one when the wave vector goes to zero. Instead, 
it approaches the canonical partition function of the model. Furthermore, note that the 
characteristic function is a Fourier Transform in a 3n-dimensional space. 

As stated before by Yamakawa|f|, the exact evaluation of the characteristic function (or 
the polymer propagator) for semiflexible chains is not possible at present. Therefore, we have 
developed a new approximation to evaluate this function. This new mathematical approach 



was developed in such a way that the most relevant physics of the problem is not altered by 
the approximation. Specifically, the proposed approach keeps the thermodynamics (partition 
function) of the model exact. Moreover, all the properties of fully flexible chains (a — > 0) and 
infinitely stiff chains (a — > oo) are preserved exactly. Consequently, this approach captures 
both asymptotic limits exactly and provides an approximate description of the crossover 
behavior. In addition, our treatment of the problem uses the exact expression of the mean- 
square end-to-end distance. Consequently, this quantity and the mean-square radius of 
gyration are exact. Another important property of our approach is that it keeps the local 
inextensibility constraint intact. Therefore, our chains have finite extensibility and this 
model can be used to compute the force-elongation relationship of semiflexible polymers. 
We describe our approximation hereafter. 

Let us start by computing the following class of integrals 

Gj = J dRjS (\Rj\ - I) exp (iRj ■ k + aR j+1 ■ R,-), (8) 

which are present in the characteristic function. The wave vector k is constant and can be 
chosen in the direction of the versor z. Writing all the vectors in spherical coordinates we 
can express Gj as follows 

Gj= Id (\Rj\) |R/ 5 (\Rj\ - I) f d6j sin (Oj) exp {i |R,-| • |k| cos (Oj)} 

J Jo (9 ) 

xexp{a|Rj| • |Rj+i| cos (Oj+i) cos (Oj)} / d(fj exp cos ((pj) + (3 sin ((fj)) , 
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where 7 and (3 are defined as follows 
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The ^-integrals can be done exactly. The result is 
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Fj (7, 0) = / d(fij exp (7 cos Upj) + f3 sin ((pj)) = 2irl (a |R,- 1 |R,-+i | |sin (Oj) sin (0 j+ i) |) , 

where Iq (x) is the Bessel function of second class[27]. After we integrate the delta function, 
the function Gj becomes 
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(12) 



We now replace this expression into Eq. (0) and obtain the characteristic function 



it n—1 

K (k, a, n, I) = 2~ n f d9 n sin (9 n ) exp {il ■ |k| cos (9 n )} f ]J d6j sin (0, 

JO JO j=\ 



(13) 



X 



exp | ^/ • |k| + al 2 cos (0j_ 



cos ( 



,)}/o(«/ 2 



| sin {Vj) sm (c^+i 



The evaluation of K (k,a,n,l) is done by iterations. First, we take the term j = 1, 
redefine |k| — > k/l and a — > a/l 2 , and remove the factor 2tt1 2 from the definition of Gj (6) 
in Eq. (jT^j) ; then, we can write 
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which is exactly doable The result is 
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The next step in the iterative process is the evaluation of G 2 (0 3 ) given by 



G 2 (0 3 ) = / rf0 2 sin(0 2 )exp{^ + «cos(03)]cos(0 2 )}/o(«sin(0 2 )sin(0 3 ))G 1 (0 2 ). (16) 
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This integral is not exactly doable. Consequently, we proceed to approximate it such that the 
asymptotic limits of flexible and stiff polymers are captured exactly. Thus, the expression 
that we obtain will give an approximate crossover behavior between the aforementioned 
limiting regimes. Note that in the limit of very stiff chains, a — > oo, all the segments will 
be parallel to each other. In other words, when a — > oo, 2 — 3 . Then, in this limit we can 
say that 



G 2 (0 3 ) - [Gi (0 3 )] 2 . (17) 

In the other limit, a — > 0, G\ (9 2 ) is independent of 9 2 . Therefore, Eq. (|TTJ) is also valid in 
the limit of flexible chains. Thus, we conclude that Eq. (|T7jl is a good approximation for 



Gi (9 3 ) since it captures the asymptotic limits exactly and provides an approximate crossover 
behavior for G2 (6*3). 

The iteration of the aforementioned approximation n — 1 times leads to the following 
expression for the characteristic function 
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Note that this expression gives the exact canonical partition function of the model 
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Thus, this first part of the approximation preserves both asymptotic behaviors and the 
thermodynamics of the problem intact. 

Let us now proceed to evaluate the approximate expression of the characteristic function, 
K (k,a,n). The integral in Eq. (fTHl) is not exactly doable thus, we evaluate it using a 
variational procedure. Let us introduce the following anzats 
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where the parameters g a , n and v a , n are determined from the constraints imposed by the 
physics of the problem as described below. One of the requirements is that the flexible and 
rigid limits are captured exactly by the model. This requires that the parameters must 
behave in the following way 
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Using Eq. (f20"|) . we can approximate the characteristic function as follows 
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Note that the term 



ikg a , n cos {0 n ) I \Ja 2 - k 2 vl nl (23) 

goes to zero as a" 1 in the limit of a — > oo and, when a — > 0, it also approaches zero because 
g a ,n ~> 0. Thus, neglecting this term does not alter the predictions of the model for the 
flexible and stiff limits. Consequently, we approximate Eq. ()2*2*1) as follows 
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exp {ik cos {0 n ) [1 + (n - 1) flf a , n ]}, 
which is exactly doable. The final expression for the characteristic function is 
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We note that the expression of the characteristic function given by Eq. (|2o|) recovers the 
exact expression of the canonical partition function of the model, Eq. 1)19)). in the limit of 
k -> 0. 

Let us now proceed to determine the values of the parameters u a>n and from the 
physics of the problem. We first look at the force-elongation behavior predicted by this 
model. This curve is given by the following mathematical expression 
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where F is the applied force and L is the average end-to-end distance of the polymer chain 
in the direction of the force. The physics of the problem imposes the following constraint 
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(27) 



which represents the finite extensibility of the polymer chain. In other words, the polymer 
chain cannot be stretched more than its total contour length. This constraint, as expressed 
by Eq. ()27)) . results in the following relationship between the parameters u a>n and g a>n 



9a,n — 1 — v a,ni (28) 

which is in perfect agreement with the required asymptotic behaviors given by Eq. f!21[) . 

Equation (|28j) gives one of the two equations required to determine the parameters v a >n 
and g a ^ n completely. The second equation is obtained from the mean-square end-to-end 
distance, (R 2 ) an - We require that our approximation reproduce this statistical quantity 
exactly. The exact mathematical expression of this average is lH l^. 
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where C (a) is the Langevin functi 

In order to derive the second relationship between v an and g a , n , we divide the character- 
istic function, Eq. 1)2 5 J) , by the canonical partition function and expand this ratio in powers 
of the wave vector k to second order. The result is the following 
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which completes our approximation. Equations (J2EJ), (jZH]) and 1)371)1 determine 2/ ajn and 
g a<n completely. Furthermore, the use of the exact expression for (R 2 ) an assures that our 



approximation predicts not only (R ) an exactly, but also \R g ) since they are related by 
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The final expression for v an is 
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and our approximation is complete. 



B. Evaluation of the Polymer Propagator using the GBT 



Replacing the expression given by Eq. (|25|) into Eq. © we obtain the following approx- 
imate expression for the polymer propagator 
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In the limit of a — > 0, Eq. (|33l) becomes 
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which is the exact expression for the polymer propagator of the Random Flight Modelj^. 
Similarly, in the limit of a — > oo, we can perform the following expansion valid for large 
values of a 

sinh{ y/a 2 - ^U}]" 1 
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and compute the polymer propagator. The result is 
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which is the polymer propagator of an infinitely stiff polymer chain |3Jj. This propagator is 
not normalized. 

We now combine the approximate expression of the characteristic function, Eq. (|25|) . 
with the Generalized Borel Transform to compute the polymer propagator and the single 
chain structure factor of the model. 

We start the evaluation of the polymer propagator by rewriting Eq. (j3~3l) as follows 

, J (R - 1 - (n - 1) g ai „, n, a) - J (R + 1 + (n - 1) g a>n , n, a) 

4vr R [1 + (n - 1) g^ n \ 
where J (r, n, a) is defined by the mathematical expression 



This integral is evaluated exactly using GBT|2H Q, Q- 

A brief summary of the GBT 

technique can be found in Appendix A. 

We present the most important steps of the calculation hereafter and leave all the math- 
ematical details for Appendix B. We first define an auxiliary function G (b, n, a) as follows 
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from which the function J (r, n, a) can be computed as the analytic continuation to the 
complex plane 



J (r, n, a) = Re {G {b = —ir, n, a)} . (39) 

Consequently, J(r,n,a) can be evaluated from the Laplace Transform given by Eq. f3~Kl) . 
Following the technique we can write 
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The integrations are computed using well-established properties and expansions of the 
functions sin (x) and [1 — x 2 ] M [28]. After some straightforward algebra, we can write the 
analytical solution of G (b, n, a) for any even number of segments as follows 
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where 3 F 2 ([,,],[,], x) and V (x) are the Generalized Hypergeometric[3J| and GammaJ 
functions, respectively . 

Replacing Eq. pTf into Eq. (j39|l and computing the analytic continuation to the complex 
plane through the substitution b = —ir we obtain 



The imaginary part of the Generalized Hypergeometric function is calculated using its well- 
known analytical properties j^. This function is an analytic function for values of the 
modulus of the argument z less than one and its continuation to the rest of complex plane 
generates a cut on the positive real axis starting at Re (z) = 1. This implies that only values 
of the argument, K n ~ 2fc ~ 2 1 ) v "."] ; larger or equal to one will contribute to the imaginary part 
of 3F2 (z). Consequently, this condition reduces the number of terms in the k— sum such 

. The explicit evaluation of Im { 3 F 2 (z)} can 
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that the last term of Eq.(JHJ is k 
be found in Appendix C. 

Finally, we replace Eq. (|C15J) into Eq. (jl2]l to obtain the exact expression for J(r, n, a) 
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Note that Eq. is just a sum of polynomials in r. The sum over the index k is the 

one obtained for the Random Flight Modelj^l and imposes the finite extensibility of the 
polymer chain. The sums over the indexes (3 and L are a consequence of the stiffness of the 
polymer backbone. The expression of J (r, n, a) can be rewritten as follows 
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which can be used for further approximation if so desired. 

The expression for J (r, n, a) given by Eq. (EH^) was derived for even number of segments 
but, its validity for odd number of segments larger than two can be proved by analytic 



continuation. Finally, in order to obtain the polymer propagator for semiflexible chains we 
have to replace Eq. (JIHJ) into Eq. (|36|) . Observe that, after the replacement, Eq. (|36|) 
recovers the known exact solution of the Random Flight Model 0] when the limit a — > is 
taken. Indeed, the only term different from zero is the one for which (3 — 0. 



Finally, we conclude our calculations of the WCM 
factor which is defined by the following formula 



3V computing the single chain structure 
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III. RESULTS AND DISCUSSION 

Figures 1, 2 and 3 show the prediction of the normalized polymer propagator, Eq. p6|). as 
function of the end-to-end distance for chains with 5, 10 and 30 Kuhn segments, and different 
values of the semiflexibility parameter a. The numerical evaluation of the propagator was 
done using Eq. (|43l) . The sum over the index (3 converges quickly even for large values of the 
semiflexibility parameter a. Indeed, even with less than 80 terms in the sum we obtained 
a relative precision of 10~ 4 . The figures clearly show that the location of the peak in the 
polymer propagator (multiplied by R 2 ) moves toward larger values of R when the stiffness 
of the polymer backbone increases. This behavior is in good qualitative agreement with 
previous results arising from computer simulation studies and theoretical approaches 
based on the Maximum Entropy Principle [20]. This is the correct result because the stiffer 
the polymer backbone, the higher the energetic penalty to bend the chain. Consequently, 
those configurations of the macromolecule with small end-to-end distance will be more and 
more hindered as the stiffness increases while those configurations with large end-to-end 
distance should be more and more favored. Therefore, the peak should shift toward larger 
values of R when the stiffness increases. 

Figure 4 shows the polymer propagator for polymer chains with 5, 10 and 30 Kuhn 
segments and a fixed value of the semiflexibility parameter a(= 3.0). This figure shows that 
the longer the polymer is, the more it behaves like a flexible chain since the location of peak 
(=end-to-end distance divided by the contour length) moves toward smaller values. In other 
words, the longer the polymer is, the less relevant the stiffness of the backbone becomes. 



Figures 5, 6 and 7 show the behavior of S (k, n, a), Eq. (j3HJ), as a function of k for different 
values of a and three values of n (5, 10 and 30). The figures clearly show a decrease of the 
single chain structure factor with increasing k until it reaches a plateau at infinite k. Note 
that our computations predict that the decrease of the single chain structure factor for small 
values of k should be faster in the case of stiff polymers than in the case of flexible ones. 
This is a consequence of the fact that rigid polymers have a larger radius of gyration than 
flexible ones for a fixed chain length. In addition, the decrease of the structure factor for 
large values of k is faster for flexible polymers than for stiff ones. Indeed, our computations 
predict that the structure factor goes as k~ e for large values of k just before the plateau 
is reached. For values of n=5, 10 and 30, the values of 6 that we got were 1, 1.08 and 
1.3 for a=10 (rigid) and 1.64, 1.9 and 2 for a=0.33 (flexible). These results are in good 
agreement with the fact that the structure factor of polymers with large chain length should 
scale as k~ d for large k values where d is the fractal dimension of the object (2 for a flexible 
polymer and 1 for a rigid polymer). Consequently, these results imply that for short polymer 
chains (n = 5), a value of a = 10 is high enough to make this polymer behave like stiff rod, 
k~ l . On the other hand, as the chain length increases we observe that a = 10 is not high 
enough to make the polymer behave as a rigid rod and deviations from the power law k~ l 
are observed. In the case of a = 0.33 our calculations predict that a value of n equal to 
five is not high enough to recover the scaling behavior of flexible chains, k~ 2 . But, as the 
number of segments increases, the exponent approaches the value of two, indicating that 
the polymer chain behaves more and more like a flexible one. This result was also showed 
in Figure 4. The figures also show that the approximation we developed in the previous 
section gives a smooth crossover behavior from the rigid to the flexible limit. 

In order to make the presentation of our work more balanced and objective, we proceed to 
compare our results with the predictions of two other models. We start with our prediction 
for the single chain structure factor, S (fc), and compare it with the expression obtained by 
Kholodenko |l7| . It has been shown that Kholodenko's result can describe experimental data 
quantitatively |l9j. Thus, a comparison between our expression for the single chain structure 
factor and Kholodenko's will help us gauge the quality of the approximations used in our 
treatment of the WCM. Figures 8 and 9 show the comparison for polymers with n = 30 
and two different values of the semiflexibility parameter (shown in the plots). We have 



checked that other values of the parameters a or a and n give quantitative agreements of 
similar quality. In all the cases studied we found that the relationship a = 2a always gives 
excellent quantitative agreement between the predictions of both models. Thus, our single 
chain structure factor should agree very well with the experimental data of Ballauff and 
coworkers. Still, since both models have different origins, some very small differences can be 
observed in the case of flexible chains (Figure 8). 

We now proceed to compare our expression of the polymer propagator with the one com- 
puted by Wilhelm and Frev[ll|. Figure 10 shows this comparison for the case of polymers 
with n — 5. The continuous curves are the results of our calculation with a = 1,3,5 and 
10. The dashed curves were constructed based on the work of Wilhelm and Frey where we 
adjusted the bending modulus such that the location of the peak in the propagator matched 
our results. This gives a better picture of the differences and similarities between both re- 
sults. Figure 10 shows quantitative agreement between both results when the stiffness is 
low. As the stiffness increases Fig. 10 shows that the qualitative behavior of both prop- 
agators is still the same. For example, both results predict that the location of the peak 
moves toward larger values of the end-to-end distance and that the distribution becomes 
narrower. The main difference between both results is quantitative in nature. Our distribu- 
tion becomes narrower than the one predicted by Wilhelm and Frey's work by a factor of two 
approximately which, in turn, generates a higher peak (the distributions are normalized). 

Finally, let us conclude this section by rationalizing the origin of the discrepancy between 
the propagators and the reason why this does not affect the structure factor significantly. We 
first note that differences in the propagators are to be expected because both, ours and Frey's, 
calculations are based on different approximations and versions of the model (continuum 
or discrete). Let us now proceed to rationalize the origin of the discrepancy between both 
predictions. Equations © and (|l5j) define the polymer propagator and single chain structure 
factor in terms of the characteristic function. Observe that the propagator is a Fourier 
transform of the characteristic function. Consequently, the oscillatory nature of the complex 
exponential generates partial cancellation of the contributions to the integral arising from 
different parts of the interval of integration. This cancellation magnifies any inaccuracies 
made in the approximation of the characteristic function. Moreover, the consequences of 
this cancellation become more and more important as the polymer becomes stiffer because 



the characteristic function itself adopts an oscillatory behavior which, in principle, is out 
of phase with respect to the complex exponential. Therefore, the stiffer the polymer is, 
the more important the consequences of the approximation become. On the other hand, 
the evaluation of the structure factor from the characteristic function does not involve any 
oscillatory function. Consequently, any small inaccuracy made in the approximation of the 
characteristic function will remain small in the expression of the structure factor, as shown 
previously. 

IV. CONCLUSIONS 

The results obtained in this paper show that the Generalized Borel Transform is a very 
useful computational tool for the statistical mechanics of single semiflexible polymer chains. 
Indeed, the results presented in this paper clearly show that GBT is able to compute polymer 
propagators for single chain problems exactly. This capability of the technique is a direct 
consequence of its mathematical simplicity (the GBT requires elements of basic calculus 
and some fundamental knowledge of complex variables). Consequently, it does not add any 
mathematical complexity to the physics of the starting model. 

Our analysis of the Wormlike Chain Model was based on an approximate expression of 
the characteristic function. The exact evaluation of this function is not possible at present. 
Therefore, we developed a new approximation that preserves the most relevant physical 
characteristics of the model intact. Specifically, our approach keeps the thermodynamics 
of the model, the flexible and rigid limits, the mean square end-to-end distance and the 
finite extensibility of the model intact while providing an approximate expression for the 
characteristic function for intermediate values of the stiffness of the polymer chain. The 
polymer propagator was obtained exactly from the approximate characteristic function using 
GBT. Note that the propagator is approximate not because of the GBT, which computes 
this quantity exactly, but because of the approximate nature of the characteristic function. 
Our expression for the propagator shows a peak that shifts toward larger values of the 
end-to-end distance as the stiffness of the polymer backbone is increased, in agreement 
with other theoretical and computational treatments of the model. We also found that, 
in the low wave vector region, the structure factor decreases faster with increasing wave 



vector when the stiffness increases. This was rationalized in terms of the behavior of the 
radius of gyration. Similarly, we found that, in the large wave vector region, the structure 
factor of flexible chains decreases faster than the one of rigid polymers. This was compared 
with the behavior of very long polymer chains whose behavior for large wave vectors is 
known exactly. We also compared of the predictions of our calculation with established 
results for the single chain structure factor and polymer propagator. Excellent quantitative 
agreement was observed between our prediction for the single chain structure factor and 
the one predicted by Kholodenko's model. The polymer propagator was compared with the 
prediction of Wilhelm and Frey. Very good quantitative agreement was observed between 
the predictions of both models for low values of the stiffness. For stiff polymers, quantitative 
deviations were observed and the origin of the deviations was rationalized. 

The proposed approach to semiflexible polymers can also address semiflexible polymers 
with other topologies like ring and m-arm star polymers. The procedure should be similar 
to the one presented in this paper but, the characteristic function will have a different 
mathematical expression. 

We conclude this section with a discussion of the characteristic function of the WCM 
which limits our ability to solve this model exactly. Let us rewrite this function. The 
expression is 

/n I n n—1 \ 

d {R fc } ft <f (|Rj| - exp » E R ^ ' k i + a E R *+i ' R * > ( 46 ) 

3=1 \ j=l i=l ) 

where we have replaced the wave vector k by a group of wave vectors kj. This expression 
is the one of the characteristic function when all the kj are equal to k. Observe that if we 
replace kj by ihj where % is the imaginary unit, then 

/n In n—1 \ 

d {R fe } J] S ( \Rj I - I) exp - ]T Rj • bj + a £ R i+1 • R , (47) 
j=l \ j=l z=l / 

which has the form of a Laplace Transform in 3n dimensions. But, the GBT computes 
Laplace Transforms very accurately or even exactly. Consequently, a generalization of GBT 
to many dimensions might lead to an exact or very accurate expression of the characteristic 
function of the WCM. This expression can be further used to compute the polymer propa- 
gator using GBT. Thus, we speculate that such extension of the GBT technique might allow 



us to solve very accurately or even exactly the Wormlike Chain Model. In general, such 
extension of GBT might allow us to solve other models of polymer chains of the form 



K ({kj} , n, parameters) = J d {R&} exp Rj • kj + H [{R;} , parameters]^ , (48) 

where H [{R;} , ...] is the Hamiltonian of the model. Thus, helical wormlike polymers and 
other models might be mathematically tractable with this generalization of the GBT. 
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Appendix A: THE GENERALIZED BOREL TRANSFORM (GBT) 

Let us briefly present the mathematical aspects of the GBT in connection with the com- 
putation of the Laplace-Mellin transform [21I Q, I3| • We start with the expression 



POO 

S(g,a } n) = / x n H (x, a) exp (— gx) dx, g>0 
Jo 

where we have explicitly extracted a factor x n from the function to be transformed. 
Defining the Generalized Borel Transform (GBT) of 5" as 



(Al) 



B x (s,a,n) 



exp [s/77] 



1 

A?7 



+ 1 



-As 



S (g,a,n)d(l/rj) , Re (s) < 



(A2) 



where A is any real positive non-zero value and I/77 = A (exp (g/X) — 1), we can invert Eq. 
(|A2J) in the following way 



OO POO 



S (g, a, n) = 2A 2 (1 — exp (— g/X)) I I exp [G (w,t, g, A, a,n)] dwdt, (A3) 

J —00 J — 00 

where the explicit expression of G (w, t, g, A, a, n) is not important for our present purposes 
(for more details see Ref. 



The expression given by Eq. (IA3|) is valid for any non-zero, real and positive value of the 

parameter A. But, the resulting expression for S (g, a, n) does not depend on A explicitly. 

Thus, we can choose the value of this parameter in such a way that it allows us to solve Eq. 
A3I). T he dominant contribution to the double integral is obtained using steepest descent 
3fil . I37T ] in the combined variables [t, w). In doing so, one first computes the saddle point 

t (g,a,n) and w Q (g,a,n) in the limit A ^> 1 and then checks the positivity condition j^Sl 

(the Hessian of G at this point should be positive) obtaining 

to = In —, — j — ■ ~t r =t (g,a,n) , w = In [x (g, a, n)} = w (g, a, n) , (A4) 

f(Xo{g,a,n),a,n)\ 

where x (g, a, n) is the real and positive solution of the implicit equation coming from the 
extremes of the function G in the asymptotic limit in A. Therefore, one obtains the following 
equation 



x o9 — f ( x o, a, n) [f (x , a, n) + 1] 



(A5) 



where 



f(x ,a,n) = 1 + n + x, 



din [H (x , a)] 

dXn 



(A6) 



In the range of the parameters where / (x , a, n) 3> 1 which is fulfilled when n>l, and 
assuming that there is only one saddle point, we can retain the first order in the expansion 
of G around the saddle point. 

Finally, we obtain the approximate expression for the starting function S (g, a, n) 



S Ap rox (g, a, n) = V2ne 1/2 



f[xo,a,n] + l r m+1 



D [x , a, n] 



\Xo 



H [x , a] exp [— / [x , a, n]) , (A7) 



where 



D (x , a, n) = -x — a ' — [1/2 + / (x a , a, n)} + / (x , a, n) [1 + / (x a , a, n)} . (A8) 

UltA/ Q 

Note that the expression given by Eq. (|A7I) is valid for functions H (x, a) that fulfill the 
following general conditions: 

1) the relation given by Eq. (|A5j) must be biunivocal. 



2) D (x , a, n) has to be positive and x c ^°' a ' n ^ — 2f (x a , a, n) has to be negative in x . 

3) / (x , a, n) 1. In particular, this condition is fulfilled when n 1. 

These conditions provide the range of values of the parameters where the approximate 
solution, Eq. (IA7jl . is valid. 

In summary, the GBT provides an approximate solution ,Eq. QA7J) . to amplitudes with 
the mathematical form given by Eq. IjAljl . The calculation consists of solving the implicit 
equation Eq. (|A5J) for n 3> 1 to obtain the saddle point and replace it into Eq. (|A7Jl . 

Let us now focus on amplitudes with the mathematical form of a Laplace transform 

/•CO 

S (g,a) = H (x, a) exp (—gx) dx g > 0. (A9) 

This kind of amplitudes can be mapped onto expressions of the form given by Eq. (I A II) . In 
order to use the GBT on Eq. (fA"9l) . we use the following relationship between Eq. (|A9|) and 
Eq. (JSH) 

S(g,a,n) = n n —S(g,a), (A10) 

which can be inverted to give 

„ „ n— 1 

S (g, a) = (-)" / dg---J dg S (g,a,n) + J2 c p ( a > n ) 9 P - (All) 



p=0 



The finite sum comes from the indefinite integrations. Note that all the coefficients vanish 
whenever the Laplace transform, Eq. (jA9j) . fulfills the following asymptotic behavior 



lim S(g, a) =0. (A12) 



In addition, the expression, given by Eq. (jAllj) . is valid for any value of n, in particular for 
n ^> 1. Consequently, if Eq. (|A12|) . is fulfilled, then the analytical solution reads 



lim r r 

S(g,a)= (-) dg--- dgS Aprox (g,a,n), (A13) 

n — > oo 



where, for n ^> 1, we can use the expression given by Eq. (|A7|) . It is important to note that 
it is the limit n — >• oo that makes the saddle point solution, Eq. (|A7J) . an exact solution for 



Eq. (IA3I) . Thus, as long as the n indefinite integrals can be done without approximations, 
as it is in our case, the result for S is exact. 



Appendix B: EVALUATION OF THE POLYMER PROPAGATOR 

We start the evaluation of the polymer propagator by rewriting Eq. (J3*%J) as follows 



G (b, n, a) 







•n— 1 



dc™- 1 



die exp (— wb) exp 

1-1 



sinh | ^Ja 2 — ifV 2 n j 



a 2 — w 2 v 2 n 



c=0' 



where GA (b, a, c) is 



c=0 



(Bl) 



GA(b,a,c)= / dw[exp(—wb)H(w,a,c)], 
Jo 

and if (w, a, c) is given by the expression 



(B2) 



H (w, a, c) = exp c 



sinh | sja 2 - w 2 ^n\ 



a 2 - w 2 z/ 2 n 



(B3) 



The integral expressed by Eq.(JB2j satisfies all the requirements of the GBT technique. 
Then, we evaluate it in the following way 



GA (b, a, c) 



lim 

iV -> oo 



db / dbGA 



(B4) 



N 



where 



GAn (6, a, c) = / dtu 
Jo 



w exp (—wb) H (w, a, c) 



(B5) 



In the asymptotic limit of iV — ► oo the GBT provides an analytical solution for Eq. (|Bf 
Following the technique, we solve the implicit equation, Eq. (jA5|) . for the saddle point w a . 
The asymptotic solution is 



+ 3/2 



iV> 1. 



(B6) 



Replacing this expression for w Q in the expression provided by the GBT, Eq. (A7), we obtain 
GA N (b, a,c)^ V + l) H (1±1H^ a} c \ iv»l. (B7) 



Furthermore, we replace Eq. (jB7jl into Eq. (|B4|) and the resulting expression into Eq. ljBlJ) . 
and we exchange the order of the operators in the resulting expression. In other words, we 
first evaluate the n derivatives with respect to c and, afterward, we take the limit of c — > to 
obtain the expression given by Eq. pTIjl . Next, we solve the N integrations using standard 
properties and expansions of the functions sin (x) and [1 — x 2 ] M 
as follows 



28l |. and write G(b,n,a) 



n — 2 
1 ^ > - 



z fc=0 



n — 1 
k 



M(N,n,k,a,b), (B8) 



where n is even and 



N — > OO r=0 r ' 

f>^~>) (-^Xv{N + l)Jdb---Jdb 



(B9) 



x 

0=0 



where (x) are the Laguerre polynomials 

Note that the only powers on b in Eq. (|B9|) that fulfill the asymptotic behavior of the 
function G (b, n, a) are those for which the condition r > (n + 2(3 — 1) is satisfied. Therefore, 
the N indefinite integrations are exactly doable. The result is 

r 1 _ T(2 + r- n-2(3) (-) N 

J db b 2+N-n+r-2p ~ f (2 + N - 71 + T - 2(3) fe2+r-n-2/3 " ( 10 ) 

Replacing Eq. t|B 1 0|) into Eq. (|B9|) and after the change of variables r = x + 2(3 + n — 1, we 
can write 



M (n, fc, a, 6) 



J oo 

/m- ^ (i (n 

" 13=0 



a oo 

2k _ 1))2/ 3+n-l (_ a 2)^^ 



i (n — 2k — l)v an 



x=0 



x- 



r(x + l)L|(0) lim AT^r(A^+l) 
T (x + 2/5 + ra) jv _, oo T(N + x + 1)' 



(Bll) 



Using the asymptotic properties of the Gamma mnction[2 



lim N X T (N + 1) 



iV^oo T(iV + x + l) 



Eq. (|B11)1 finally reads 



(B12) 



M (n, k, a, b) = -Im J2J2(i(n-2k- 1)) 



2/9- 



n _i (-a 2 )' 3 fi(n-2k -1) v a . 



/3=0 x=0 



r (x + 1) r (2±a + /3) 
r (x + 2/3 + n) r (2±?) 

The sum over x is exactly doable. The result is 



(B13) 



M (n, fc, a, 6) = -im £ (i (ra - 2fc - 1)) 



2/3+n-l 



—a 



2\P 



0\ 



-FD (p,n,k,b) 



where 



(B14) 



FD {P,n,k,b) 



r(/3 + l) 3 F 2 ([l A /3 + ll,[/3+2±!,/3+2 



[(n-2fc-l)i/ a ,n]^ 
b 2 



T (2/5 + n) 

2z (n - 2* - 1) T (/? + |) 3^2 ([/? + |, 1, l] , [/? + f + 1, /5 + § + |1 - ^^-^ 



T {2(3 + n + 1) 



(B15) 



Equation ()B14f) clearly shows that the imaginary part affects only the real part of the 
function FD (/3, n, k, b). Thus, the final expression for G (b, n, a) is given by Eq. (JUJ. 



Appendix C: EVALUATION OF THE IMAGINARY PART OF THE GENERAL- 
IZED HYPERGEOMETRIC FUNCTION 



We evaluate Im { 3 F 2 (z)} using the following integral representation of the Hypergeomet- 
ric function! 



3 r 2 



n+1 _ n 



in -2k- 1) 



&2(/3+l) 



6 2 



(n - 2Jfc - 1) 



2/3+n-l 



n-2fc-l 



S(l,2/3 + n- r 



[n — 2fc — 1 — x 



2/3+n-2 



x 2 + b 2 



(CI) 



dx, 



which is valid for values of n equal or larger than two. The analytic continuation to the 
complex plane is done as before through the replacement b = —ir. Then, the imaginary 
part of the Hypergeometric function is 



Im { 3 F 2 



1,2^ + ! 



' n + 1 . n 

r 2(/3+l) 



(n-2k- T 



>2/3+n-l 



B (1,2/3 + n- l)(n-2k- 1 

xlm [ [n — 2k — 1 — x] 2l3+n ~ 2 f 

Jo 



(C2) 



2 2 

x — r 



dx. 



Thus, we have to evaluate the expression 



L = I in I [n — 2k — 1 — x 



213+71 2 \ x — r] 13 1 \x + r] 13 1 



dx. 



(C3) 



After analyzing the analytical behavior of the integrand we concluded that we can ex- 
change the operations of integration and imaginary part to obtain 



L 



n-2k-\ 



[n — 2k — 1 — x] 2/3+n 2 [x + r] 13 1 /m|[x — r] 13 x | dx. 



(C4) 



Thus, we have to compute 



LS = Im 



1 



(C5) 



n. 



^ [x — r) J 

first and, afterward, we have to solve the integral given by Eq. (|C4jl . 

The analytical behavior of the function (x — i?) is well known |28l|. It is an analytic 
function for \x\ > R but, its analytic continuation to the complex plane generates a cut on 
the real axis in the range —R < Re (x) < R which provides its imaginary part. Writing 



{x - R) 



[X 



(C6) 



for integer values of (3 and x > R, and using the following integral representation 



1 

y/x — R -'a 

valid for x > R, we obtain the result 



dyexp —y\/x — R 



(C7) 



Im 



1 

(x-R) 



13+1 



(R-x) 



Re \ du exp [— iu (R — x)\ 



(C8) 



where x < R. 

Using parity's properties of the function cos (9) , we extend the range of the integration 
in Eq. (JOE!) to obtain 



Im 



V& (x-R) = — W-t- =p6 (x-R). (C9) 



{x-Rf +1 \ {x-Rf"^ " (31 d(x-R) 
Finally, we replace the expression given by Eq. (|C9() into Eq. ()C4|) and perform the 
change of variables y = x — r to obtain 



(3\ 

where Fpk (y, n, r) is defined as 



k(-) P fn-2k-l-r 

L = — 7^— I F/3 k (y, n, r) -—5 (y) dy, 



(CIO) 



F pk (y, n, r) = [n - 2k - 1 - r - yf p+n 2 [y + 2r] p 1 . 
Integrating by parts (3 times, we obtain the following final expression 



(Cll) 



L = — 



7r_ / d p F 0k (y, n, r) 

a 



P\ 1 dy 
The f3 derivatives are computed as follows 



(C12) 



d^F^ (y, n, r) 
dy? 



■ =E 

y=0 L=0 



a 1 



dyP- 



1 dy L 

[y + 2r]- p - 1 



[n — 2k — r — 1 — y] 



2/3+n-2 



y=o 



(C13) 



to obtain 



d^Fpk (y, n, r) 
dyP 





> =E 

y=0 L=0 



(3 \ (-f (2(3 - L)\ {2(3 + n - 2)\ [n - 2k - r - 1] 



23+n-2-L 



L 



{2(3 + n-2- L)\{(3)\{2r) 



2/3+l-L 



(C14) 

Finally, we replace Eq. (|C14|) into Eq. (IC12|) and the resulting expression into Eq. (IC2I) 
to obtain the final expression of Im {3F2} 



Im{ 3 F 2 } 



(3 

l)2/3+ n-l [r {(5 E 



vrr (2/3 + n) 



(n-2A; 
r ^[ n _2A;-r-l] 2/3+ "- 2 - L 
(2/3 + n-2-L)! 




(2/3 -L)! 

(2)2/3+ l-L 



(C15) 



n > 2. 



[1] M. G. Poirier et al., Phys. Rev. Lett. 86, 360 (2001); J. F. Leger et at., Phys. Rev. Lett. 83, 
1066 (1999); V. Parpura and J. M. Fernandez, Biophys. J. 71, 2356 (1996); M. D. Wang et 
al, Biophys. J. 72, 1335 (1997); A. D. Mehta, K. A. Pullen and A. Spudich, FEBS Lett. 430, 
23 (1998); A. D. Mehta, M. Rief and J. A. Spudich, J. Biol. Chem. 274, 14517 (1999),T. E. 
Fisher, P. E. Marszalek and J. M. Fernandez, Nat. Struct. Biol. 7, 719 (2000); T. E. Fisher et 
al., Trends Biochem. Sci. 24, 379 (1999); M. Carrion- Vazquez et al., Prog. Biophys. Mol. Biol. 
74, 63 (2000); T. T. Perkins, S. R. Quake, D. E. Smith and S. Chu, Science 264, 822 (1994); 
S. B. Smith, Y. Cui and C. Bustamante, Science 271, 795 (1996); T. T. Perkins, D. E. Smith 
and S. Chu, Science 276, 2016 (1997); A. F. Oberhauser, P. E. Marszalek, H. P. Erickson and 
J. M. Fernandez, Nature (London) 393, 181 (1998); F. Oesterhelt, M. Rief and H. E. Gaub, 
New J. Phys. 1, 6.1 (1999); H. B. Li, W. K. Zhang, W. Q. Xu and X. Zhang, Macromolecules 
33, 465 (2000). 

[2] J. F. Marko and E. D. Siggia, Macromolecules 28, 8759 (1995); K. Kroy and E. Frey, Phys. 

Rev. Lett. 77, 306 (1996). 
[3] G. Porod, Monatsh. Chem. 80, 251 (1949); O. Kratky and G. Porod, Rec. Trav. Chim. 68, 

1106 (1949). 

[4] L. Eichinger et al, Biophys. J. 70, 1054 (1996); P. Janmey Cell Membranes and the Cy- 
toskeleton, Vol. 1A of Handbook of Biological Physics, chapter 17, page 805 (North Holland, 
Amsterdam, 1995); P. A. Janmey, Curr. Op. Cell. Biol. 2, 4 (1991). 

[5] N. Saito, K. Takahashi and Y. Yunoki, Journal of the Physical Society of Japan, 22, 219 (1967). 

[6] H. Yamakawa, Helical Worm-like Chains in Polymer Solutions (Springer, Berlin, 1997). 

[7] M. Fixman and J. Kovac, J. Chem. Phys. 58, 1564(1973). 

[8] R. A. Harris and J. E. Hearst, J. Chem. Phys. 44, 2595 (1966). 

[9] H. E. Daniels, Proc. Roy. Soc. Edinburgh A63, 290 (1952); W. Gobush, H. Yamakawa, W. H. 



Stockmayer and W. S. Magee, J. Chem. Phys. 57, 2839 (1972). 
[10] H. Yamakawa, Pure and Appl. Chem. 46, 135 (1976). 

[11] J. Wilhelm and E. Prey, Phys. Rev. Lett. 77, 2581 (1996); P. L. Hansen and R. Podgornik, J. 

Chem. Phys. 114, 8637 (2001). 
[12] M. G. Bawendi and K. F. Freed, J. Chem. Phys. 83, 2491 (1985); J. B. Lagowski, J. Noolandi 

and B. Nickel, J. Chem. Phys. 95, 1266 (1991); G. A. Carri and M. Muthukumar, J. Chem. 

Phys. 109, 11117 (1998). 
[13] R. G. Winkler, P. Reineker and L, Harnau, J. Chem. Phys. 101, 8119 (1994). 
[14] C. M. Marques and G. H. Fredrickson, J. Phys. II France 7, 1805 (1997). 
[15] L. Harnau, R. G. Winkler and P. Reineker, J. Chem. Phys. 104, 6355 (1996); K. Kroy and E. 

Frey, Phys. Rev. E 55, 3092 (1997). 
[16] A. L. Kholodenko, Ann. Phys. 202, 186 (1990); A. L. Kholodenko J. Chem. Soc. Faraday 

Trans. 91, 2473 (1995); A. L. Kholodenko, Phys. Lett. A 141, 351 (1989); A. L. Kholodenko, 

Phys. Lett. A 178, 180 (1993); A. Kholodenko and T. Vilgis Phys. Rev. E 50, 1257 (1994), A. 

L, Kholodenko and T. A. Vilgis Phys. Rev. E 52, 3973 (1995). 
[17] A. L. Kholodenko, Macromolecules 26, 4179 (1993). 

[18] A. Kholodenko, M. Ballauff and M. Aguero Granados, Physica A 260, 267 (1998). 
[19] M. Ballauff et al, Macromolecules 30, 273 (1997). 
[20] R. G. Winkler, J. Chem. Phys. 118, 2919 (2003). 

[21] L.N.Epele, H. Fanchiotti, C.A. Garcia Canal and M. Marucho, Phys. Lett. B523, 102 (2001). 
[22] L.N.Epele, H. Fanchiotti, C.A. Garcia Canal and M. Marucho, Phys. Lett. B556, 87 (2003). 
[23] L.N.Epele, H. Fanchiotti, C.A. Garcia Canal and M. Marucho, Nucl. Phys. B583, 454 (2000). 
[24] S. M. Bhattacharjee and M. Muthukumar, J. Chem. Phys. 86, 411(1987). 
[25] K. F. Freed in Advances in Chemical Physics, Vol. 22, Pages 1-128, I. Prigogine and Stuart A. 

Rice Editors (John Wiley and Sons, New York, 1972). 
[26] G. Arfken, Mathematical Methods for Physicists (Academic Press, Inc., New York, 1985). 
[27] M. Abramowitz and I. Stegun, Handbook of Mathematical Functions (Dover, New York, 1970). 
[28] I.S. Gradshteyn and I.M. Ryzhik, Table of Integrals, Series, and Products (Academic Press, 

New York, 2000). 
[29] H. Benoit and P. Doty, J. Phys. Chem. 57, 958 (1953). 



[30] K. F. Freed, Renormalization Group Theory of Macromolecules (John Wiley & Sons, New 
York, 1987). 

[31] P. G. de Gennes, Biopolymers 6, 715 (1968). 

[32] H. Bateman, Higher Transcendental Functions (Mc Graw Hill, New York, 1953). 

[33] M. C. Wang and E. Guth, J. Chem. Phys. 20, 1144 (1952); C. Hsiung, H. Hsiung and A. 
Gordus, J. Chem. Phys. 34, 535 (1961). 

[34] W. L. Mattice and U. W. Suter, Conformational Theory of Large Macromolecules: The Ro- 
tational Isomeric State Model in Macromolecular Systems (John Wiley & Sons, New York, 
1994). 

[35] J. Wilhelm and E. Frey, Phys. Rev. Lett. 77, 2581 (1996). 

[36] N. Bleistein and R. Handelsman, Asymptotic Expansion of Integrals (Dover, New York, 1986). 
[37] E. T. Copson, Asymptotic Expansions (Cambridge University Press, 1965). 
[38] H. Jefrey and B. S. Jeffrey, Methods of Mathematical Physics (Cambridge University Press, 
1966). 



List of Figures 



FIG. 1: Normalized polymer propagator 4nR 2 P (R,n,a) /Z n (a) versus R/n for n = 5. Con- 
tinuous line (a = 0.33) , dotted line (a = 1.0), dashed line (a = 3.0), long dashed line 
(a = 5.0) and dashed-dotted line (a = 10.0). 

FIG. 2: Normalized polymer propagator 4nR 2 P (R, n, a) /Z n (a) versus R/n for n = 10. Con- 
tinuous line (a = 0.33) , dotted line (a = 1.0), dashed line (a = 3.0), long dashed line 
(a = 5.0) and dashed-dotted line (a = 10.0). 

FIG. 3: Normalized polymer propagator 4ttR 2 P (R, n, a) jZ n (a) versus R/n for n = 30. Con- 
tinuous line (a = 0.33) , dotted line (a = 1.0), dashed line (a = 3.0), long dashed line 
(a = 5.0) and dashed-dotted line (a = 10.0). 

FIG. 4: Normalized polymer propagator 4nR 2 P (R, n, a) /Z n (a) versus R/n for a = 3.0. Con- 
tinuous line (n = 30), dashed line (n = 10) and dashed-dotted line (n = 5). 

FIG. 5: Single chain structure factor S (k, n, a) versus wave vector k for n = 5. Continuous line 
(a = 0.33), dotted line (a = 1.0), dashed line (a = 3.0), dashed-dotted line (a = 10.0) 
and circles (best fits to the power law in the range k G (2, 3) for a = 0.33 and a = 10.0). 

FIG. 6: Single chain structure factor S(k,n, a) versus wave vector k for n = 10. Continu- 
ous line (a = 0.33), dotted line (a = 1.0), dashed line (a = 3.0), dashed-dotted line 
(a = 10.0) and circles (best fits to the power law in the range k e (2, 3) for a = 0.33 
and a = 10.0). 

FIG. 7: Single chain structure factor S(k,n,a) versus wave vector k for n = 30. Continu- 
ous line (a = 0.33), dotted line (a = 1.0), dashed line (a = 3.0), dashed-dotted line 
(a = 10.0) and circles (best fits to the power law in the range k G (2, 3) for a = 0.33 
and a = 10.0). 

FIG. 8: Single chain structure factor S(k,n, a) versus wave vector k for n = 30. (Line) 
Kholodenko's model with a = 1, (points) this work with a = 2. 

FIG. 9: Single chain structure factor S(k,n,a) versus wave vector k for n = 30. (Line) 
Kholodenko's model with a = 50, (points) this work with a = 100. 



FIG. 10: Normalized polymer propagator P(R,n,a) versus end-to-end distance R in units of 
Kuhn length for n = 5. (Continuous lines) this work, (dashed lines) Wilhelm and 
Frey's results. 
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FIG. 3. Normalized polymer propagator AitR 2 P (R, n, a) /Z n (a) versus R/n for n = 30. 
Continuous line (a = 0.33), dotted line (a = 1.0), dashed line (a = 3.0), long dashed line 
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FIG. 5. Single chain structure factor S(k,n,a) versus wave vector k for n = 5. Continuous line 
(a = 0.33), dotted line (a = 1.0), dashed line (a = 3.0), dashed-dotted line (a = 10.0) and circles 
(best fits to the power law in the range k £ (2, 3) for a = 0.33 and a = 10.0). 
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FIG. 6. Single chain structure factor S (k,n,a) versus wave vector k for n = 10. Continuous line 
(a = 0.33), dotted line (a = 1.0), dashed line (a = 3.0), dashed-dotted line (a = 10.0) and circles 
(best fits to the power law in the range k 6 (2, 3) for a = 0.33 and a = 10.0). 
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FIG. 7. Single chain structure factor S (k,n,a) versus wave vector k for n = 30. Continuous line 
(a = 0.33), dotted line (a = 1.0), dashed line (a = 3.0), dashed-dotted line (a = 10.0) and circles 
(best fits to the power law in the range k 6 (2, 3) for a = 0.33 and a = 10.0). 
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FIG. 8. Single chain structure factor S (k,n,a) versus wave vector k for n = 30. (Line) 
Kholodenko's model with a = 1, (points) this work with a = 2. 
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FIG. 9. Single chain structure factor S (k,n,a) versus wave vector k for n = 30. (Line) 
Kholodenko's model with a = 50, (points) this work with a = 100. 
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FIG. 10. Normalized polymer propagator P(R,n,a) versus end-to-end distance R in units of 
Kuhn length for n = 5. (Continuous lines) this work, (dashed lines) Wilhelm and Frey's results. 
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